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1. Introduction 

At the present time about twenty Kl n~^n~e~^e~ events have been observed 
and a detailed experimental study of this decay mode will be possible in future exper- 
iments [Q. The Kl 7r+7r~e"'"e~ weak decay amplitude is dominated by the process, 
Kl — ^ tt"*" n~ 7* n~^n~ e~^e~ , where a single virtual photon creates the e~^e~ pair. 
This one photon contribution to the decay amplitude has the form 

M(i-^) = [iGe^^^y+x P-p q. + F+pI + F_p^_ ] ■u{k.Mk+) , (1.1) 

where Gp is Fermi's constant, a is the electromagnetic fine structure constant, si ^ 0.22 
is the sine of the Cabibbo angle and / ~ 132 MeV is the pion decay constant. The tt"*" and 
n~ four-momenta are denoted by p+ and p- and the e"*" and e~ four-momenta are denoted 
by /c-i- and k-. The sum of the electron and positron four-momenta is q = k- + k+. The 
Lorentz scalar form factors G, F± depend on the scalar products of the four-momenta q,p+ 
and p- . Neglecting CP nonconservation, under interchange of the pion four-momenta 

p+ P- and p- p+ (1-2) 

the form factors become 

G^G , F+^F_ , F_^F+ . (1.3) 

In this paper we compute the CP conserving contribution to the form factors G, F± us- 
ing chiral perturbation theory at one-loop order (the 0{p'^) amplitude vanishes). The 
coefficients of some of the local operators appearing at the same order in the chiral ex- 
pansion (i.e., order p'^ counter terms, where p is a typical momentum) are determined by 
the experimental value of the pion charge radius and the measured 7r~^e~^e~ and 

Kl 7r"'"7r~7 decay rates and spectra. 

We also compute (in chiral perturbation theory) an important tree level contribution 
to the form factors F± that arises from the small CP even component of the Kl state. 
This contribution to the F± form factors from indirect CP nonconservation has the oppo- 
site symmetry property under interchange of pion momenta when compared with the CP 
conserving contribution to F± (see eqs. ( pT^ ) and ( p..3|) ) . If 

— * p- and p- p^ (1-4) 
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then the CP violating one-photon form factors become 



F+ -F_ 



F_ 



(1.5) 



The decay amphtude that foUows from squaring the invariant matrix element in eq. 
(1.1) and summing over e"*" and e~ spins is symmetric under interchange of e"*" and e~ 
momenta, k- k+. Physical variables that are antisymmetric under interchange of the 
e'^ and e~ momenta arise from the interference of the short distance contributions (Z- 
penguin and 1^-box diagrams) and the two photon piece with the one photon amplitude 
given in eq. (1.1). 

In the minimal standard model the coupling of the quarks to the VF-bosons has the 

form 



Here repeated generation indices j,k are summed over 1,2,3 and g2 is the weak SU(2) 
gauge coupling. F is a 3 x 3 unitary matrix (the Cabibbo-Kobayashi-Maskawa matrix) 
that arises from diagonalization of the quark mass matrices. By redefining phases of the 
quark fields it is possible to write V in terms of four angles 61,62, 6s and 5. The 9j are 
analogous to the Euler angles and 5 is a phase that, in the minimal standard model, is 
responsible for the observed CP violation. Explicitly 



where Ci = cos 6i and Si = sin^j. It is possible to choose the Oj to lie in the first quad- 
rant. Then the quadrant of S has physical significance and cannot be chosen by a phase 
convention for the quark fields. A value of S not equal to or tt gives rise to CP violation. 

The short distance W-hox and Z-penguin Feynman diagrams depend on the Vts el- 
ement of the Cabibbo-Kobayashi-Maskawa matrix. It is very important to be able to 
determine this coupling experimentally. In this paper we calculate the contribution to 
the Kl TT~^7T~e~^e~ decay amplitude arising from the Z-penguin and W-hox diagrams 
which can be determined using chiral perturbation theory since the left handed current 
57^(1 — 75) d is related to a generator of chiral symmetry. At the present time all observed 
CP nonconservation has its origin in — mass mixing. A CP violating variable can 
be constructed in the decay Kl — > 7r"'"7r~e"''e~ that gets an important contribution from 



(1.6) 




(1.7) 
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CP nonconservation in the Z-penguin and VF-box diagrams, that is, direct CP violation. 
The variable (in the Kl rest frame) 



=< " "^^^ ■ ^^ " ^^ > , (1.8) 

\{P- xp+) ■ {k--k+)\ 

is even under charge conjugation and odd under parity. It is also odd under interchange of 
/c_(_ and The real and imaginary parts of Vts are comparable, and hence the CP con- 
serving and CP violating parts of the Z-penguin and VF-box diagrams are of roughly equal 
importance. Acp gets a significant contribution from this direct source of CP nonconser- 
vation. In this paper we calculate Acp in the minimal standard model but unfortunately 
we find that it is quite small; |Acp| ^ 10~^. 

The decay 7r+7r~e"'"e~ has been studied previously by Sehgal and Wan- 

ninger ^ and by Heiliger and Sehgal 0. These authors adopted a phenomenologi- 
cal approach, relating the Kl n~^n~e~^e~ decay amplitude to the measured Kl 
7r"'"7r~7 decay amplitude. In the systematic expansion of chiral perturbation theory we 
find important additional contributions to the Kl n~^n~e~^e~ decay amplitude for 
q'^ = {k- -f /c_|_)^ >> 4mg that were not included in this previous work. It was pointed 
out in refs. 0] and that indirect CP nonconservation from K^ — K^ mixing gives an 
important contribution to the Kl n^n~e~^e~ decay rate and consequently there is a CP 
violating observable, Bcp, that is quite large. We reexamine Bcp using the form factors 
determined in this paper. 



2. The One-Photon Amplitude 

Chiral perturbation theory provides a systematic approach to understanding the one- 
photon part of the Kl 7r+7r~e+e~ decay amplitude. It uses an effective field theory that 
incorporates the SU(3)LxSU(3)i<; chiral symmetry of QCD and an expansion in powers of 
momentum to reduce the number of operators that occur. In the chiral Lagrangian the 
tt's , -ftT's and ry are incorporated into a 3 x 3 special unitary matrix 

f2iM\ , , 

E = exp — ^ , (2.1) 



/ 



where 




M=\ TT- -7vyV2 + v/V6 . (2.2) 
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At leading order in chiral perturbation theory / ~ 132 MeV is the pion decay constant. 
Under SU(3)lX SU(3)ii; transformations the E field transforms as 

E LER^ , (2.3) 

where L e SU{3)l and R e SU{3)r. 

At leading order in chiral perturbation theory (i.e., order p^, where p is a typical four- 
momentum) the strong and electromagnetic interactions of the pseudo-Goldstone bosons 
are described by the chiral Lagrange density 
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= ^Tr{D^T:D^T:^) + vTr{mqE + rUqE^) , (2.4) 
8 

where f is a parameter with dimensions of mass to the third power and ruq is the quark 
mass matrix 


















rUs 



mg = \ md . (2.5) 



In this paper we neglect isospin violation in the quark mass matrix and set rriu = md- In 
this approximation the and have equal masses which we denote by m^, and the 
Gell-Mann-Okubo mass relation 

3m2 - Am\ + = , (2.6) 

holds. 

The effective Lagrangian for AS" = 1 weak nonleptonic decays transforms as 
{^lAr) + (27L,li?) under SU(3)L®SU(3)fl. The (8^,1^) amplitudes are much larger 
than the {21^,1^1) amplitudes and so we will neglect the (27^,17^) part of the effective 
Lagrangian. The effective Lagrangian for weak radiative kaon decay is obtained by gauging 
the effective Lagrangian for weak nonleptonic decays with respect to the U(1)q of elec- 
tromagnetism. At leading order in chiral perturbation theory the AS" = 1 transitions are 
described by 

^U) ^ asGpSif ^^ T] + h.c. . (2.7) 

4v 2 

The matrix T in ( p77|) projects out the correct flavour structure of the octet 



T= , (2. 
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and Qs is a constant determined by the measured Ks tt+tt decay rate; \gs\ — 5.1. In 
( p^) and ( p.7|) represents a covariant derivative: 

L>^E = a^E + zeA^[Q,E] , (2.9) 

where 

/2/3 \ 
Q = [ -1/3 , (2.10) 

\ -1/3 / 

is the electromagnetic charge matrix for the three hghtest quarks, u, d and s. 
The Kl state 

\Kl >^ \K2 > +e\Ki > , (2.11) 



is mostly the CP odd state 



\K, >= ^{\K' > +\k' >) , (2.12) 
with small admixture of the CP even state 

|Ki >= ^(1^0 > -|Z%) . (2.13) 

The parameter e characterizes CP nonconservation in — mixing. At leading order 
in chiral perturbation theory the tt"*" tt" 7* n~^n~e~^e~ decay amplitude 

arises though the CP even component of Kl. Writing the form factors contributing to 
Kl 7r+7r~7*asa power series in the chiral expansion 

F± = + F^^ + ... , G = G^^^ +G^^^ + ... , (2.14) 

where the superscript denotes the order of chiral perturbation theory, we find that the 
Feynman diagrams in fig. |l] give 



G'(i) = 
F 



(i)_ 32gsP{ml-ml)n\ 



+ [q^ + 2q-p+] ■ (2.15) 

^(1) _ ^ 32g8pim'j^ -rnDit^e 
[q'^ + 2q-p_] 

Despite the fact that e ~ 0.0023 e*^^ (in a phase convention where the K^ 7rn{I = 0) 
decay amplitude is real), is small it is important to keep this part of the decay amplitude. 
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Other contributions not proportional to e don't occur until higher order in chiral pertur- 
bation theory. We neglect direct sources of CP nonconservation in the one-photon part of 
the decay amplitude. Experimental information on e' suggests that they are small. 

At the next order in the chiral expansion the form factors G^'^\ arise from 
local operators and from one-loop Feynman diagrams involving vertices from the leading 
Lagrange densities in ( |2.4|) and (|2.7|) . However, the form factor G^"^^ arises solely from 
local operators as the one loop Feynman diagrams and tree graphs involving the Wess- 
Zumino term [|5| do not contribute. The contribution of the 0{p'^) local operators to 
G''^) is fixed by the measured 7r+7r~7 decay rate to be 

1^(2)1-40 . (2.16) 

The experimentally observed 7r"'"7r~7 Dalitz plot suggests that the form factor G 

has significant momentum dependence. This indicates that G*-^-* is not negligible, and that 
our extraction of G*^^-* from the rate is not completely justified [Q. 

The form-factors F^'^ get contributions both from local operators of 0{j>^) and 
from one-loop diagrams involving vertices from the leading Lagrange densities in ( |2.4| ) and 
( p.7| ) . For Kl 7r"'"7r~e"'"e~ the local operators that contribute are 

4'^ = ~'"^^-J^^ F^-rr [Q(D^ED.Et + D,T.^D^T.)] , (2.17) 



and 



+a2iij)F^"'Tr[Q{ED^E^)T{ED,j:^)] 

+a3{fi)F^"'Tr[T[Q, E]D^E^ED,E^ - TD^ED^E^Ept, Q]] 
+a^{^i)F^'^Tr[TT.D^T.^[Q,T]D^T.^] ] + h.c. 

(2.18) 

The coefficients Acr, 01,^25^3 and 04 depend on the renormalization procedure used and 



we employ dimensional regularization with MS subtraction. The dependence of the co- 
efficients Acr-,^! 2,3,4 on the subtraction point fx cancels that coming from the one-loop 
diagrams. Note that the basis of operators in eq. ( p.l8|) is slightly different than that used 
in 1^ . With this basis of operators the combination of counterterms 

WL = a3-a4 , (2.19) 
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is independent of the subtraction point ^ at one loop. 

The value of Acr is fixed by the measured 7i~^ charge radius; < >= 0.44 ± 0.02 fm^. 



The one-loop diagrams in fig. ^ give (using MS) 

XM = -(^^^f<rl> - ^^{2ln{mlli?)^ln{m\li?)\ , (2.20) 

which implies that (at the subtraction point n = IGeV) 

Acr(lG'e"l/) = -0.91 ±0.06 . (2.21) 

A linear combination of ai and a2 is fixed by the measured n^e~^e~ decay 

amplitude. Fortunately it is the same combination of ai and a2 that enters into the 
Kl — > 7r+7r~e"'"e~ decay amplitude. The one-photon part of the Tr+e+e" decay 

amplitude can be written in terms of a single form factor f{q^) 

M^^^\K+ TT+e+e-) = ^ ^ f{q') p^^ u{k.)^^v{k+) . (2.22) 



The one- loop diagrams in fig. ^ and the operators in (|2.17| ) and ( |2.18| ) give [^| 



+ + 2a2(;u)) - 4A,,(;u) + ^ 



(2.23) 



where 



= 2g8{(t)K{q^)+(t)M^) + w+ ) 
■,iq^) = dx(^-x{l-x)^ln(^- i^x{l-x)^ . (2.24) 



This relation defines the ^ independent constant [|T0[ which has been experimentally 
determined to be [ |TT 

«;+ = 0.89t^:24 . (2.25) 
Using the central values of \cr{)-GeV) and wj^ we find that 

ai{lGeV) + 2a2{lGeV) = -6.0 , (2.26) 

with an associated error around 10% (which is correlated with the uncertainty in 
\cr{^GeV)). Throughout the remainder of this work we will use the central values of 
\cr{)-GeV) and ai(lGeV) + 2a2{lGeV) and suppress the associated uncertainties. Note 
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that the contributions Xcri^GeV) and {ai+2a2){lGeV) to /(g^) are separately quite large 
but they almost cancel against each other. 

At 0{p'^) the form factors F^^ for Kl Tv~^Tv~e'^e~ decay follow from the Feynman 
diagrams in fig. ^ and tree level matrix elements of the operators in ( p.l7| ) and ( p.l8| ) . We 



find, using MS subtraction, that 



.(2) 



/ 2 2 



-4 / dx 
'o 



q^x{l — x)in 



— q^x{l — x) 



m^in I 1 



q'^x{l — x) 



mt 



+2q' 



q'^ + 2q- (p+ +p. 



1 



where 



(2.27) 



1 



2 \2 



(q^ + 2q-p-j JO 
.1 



1 ri-y I 

dyy / dx — 
^/o A^i 

(ix£n ( 1 — 



x{l — x){q^ + 2q ■ p- 



m 



— x) + m^x — m^x(l — x) 



m 



\{l-y)+ m'^y - mly{l - y) 



nl nl-y 

+3 dy dxin ( 1 

Jo Jo 



{q'^x{l — x) + 2q ■ p+xy) 



+ f dyy [ \ 
Jo Jo l^l 

+ J dx {1 + X + 



-y) + - 'mly{l - y) 
(g(l - x) +p-y) ■ (4g + 6p+ + 4p_) - 2mj^ - 2p+ ■ {q + p-) 
(x - l)(3m|. - 2m^) \ / x{l - x){q'^ + 2q ■ p-) 



+ 2p- ■ q 



£n 1 - 



m\{l — x) + m^x — m^x(l 



X 



(2.28) 



with 



IJ^i = 'mKi^-y)+'mly-mly{l-y)-q'^x{l-x)-2q-p-xy . (2.29) 
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The other function is 

c^K. = (m^ -m^){ I dy dxin - ^2^^ ^ ^2 (i _ 







mj^y + ml{l - yY 
{P+ +p-+q)- {q{l - x) + yp- - p+) 



'0 Jo 

+2 / dyy dx ^ 

Jo Jo 1^2 

/\ / m|^(x-l)\ / x(l- x)(q'^ + 2q-p. 
+ dx {1 + X+ \in\l ^ ^ ^ 







q^ + 2q-p_J \ rn^il — xY + m\-x 



(2.30) 



/i| = m^j^y + m^(l — yY — q^x{l — x) — 2q ■ p-xy . (2-31) 



with 



The GeU-Mann-Okubo mass formula (|2.6D has been used to simphfy some of the de- 

pendence on the pseudoscalar masses in (lOSD and (lOOD . F\:> is obtained from (^) 

by taking p+ p- and p- ^ p+- Notice that the combination of coefficients (ai + 2a2) 

(2) 

and Acr that appear in the expression for has a relative sign difference compared to 
the combination that appears in the expression for f{s) given in eq. (|2.23| ). The uncer- 
tainty in Xcri^GeV) and gives rise to about a 10% uncertainty in the combination of 
counterterms that appears in ( p.27|) . The one photon part of the Tv^Tv~e~^e~ decay 

amplitude is the largest and dominates the rate. In the next section we use the form fac- 
tors calculated here to obtain dr{KL — > n^'K~e~^e~) /dq'^ . One (scale independent) linear 
combination of counterterms wl = as — cla is not determined by the present experimental 
data and consequently we cannot predict the rate for — > n^n~e~^e~ . However, this 
is the only undetermined constant and the entire function dr{KL n~^n~e^e~)/dq'^ is 
experimentally accessible. 



3. The Differential Decay Rate 

The Kl n^n~e^e~ decay rate is obtained by squaring the invariant matrix element 
(1.1), summing over the e+ and e~ spins, and integrating over the phase space. Since the 
e"*" and e~ four momenta only occur in the lepton trace, Tr [^_7;y^^7^] , the phase space 
integrations over k- and /c+ produce a factor 

IiSk^_I(Sm^''^'''^'-'--'*^^'^^-''^^-''''^ . ,3.1) 

1 
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The remaining phase space integrations can be taken to be over and the sum and 
difference of the pion energies in the rest frame, Es = p\ + p^^, Ed = — p^_. The 
contribution of the form factors F± and G to dV/dq^ do not interfere. Therefore, we can 
write 

rir dVa dVp , , 

+ ^ , (3.2) 



dq' 



dq' 



where 

dTo _ 
dq^ mKp2^{'2ny^'^ 



GW4 



dEs / dEolGl 



m^Q -m^{P--(l) -m^{p+-Q) + ^P+ ■ P-){q ■ P+){q ■ P-) - Q {p+ ■ P- 



dr, 



Gla'sl 



dq^ mKp2^{2ny3q^ 



dEs / dEo 



|F+ q- p+ + F- q- p. 



q^ (iF+pm^ + iF.pm^ + 2Re{F+Fl)p+ ■ p_) 



(3.3) 



In eq. (|3.2|) and eq. (^3) difference of pion energies is integrated over the region 



E 



(max) 
D 



< Ed < E 



(max) 
D 



where 



E 



'max) 
D 



' 2mKEs + q^ -m\- Ami n 2 

— „ , 2 ^ 2 V - EsY - r 

2mKEs + q^ - mj^ 



(3.4) 



and the sum of pion energies is integrated over the region 

^{m^n) ^ < 4"""^ ^^ere 



the boundaries are 



T-\imax) I 9 

E^g ' = rriK — V 9 



E 



(min) _ mj^ -q^ +Am 



(3.5) 



2mK 



The scalar products appearing in the expression for the rates are easily expressed in terms 
of Es, Ed and q^: 



P+-P- = ^{q"^ - 2to^ + 2mKEs) 



1 

2' 
1 



q-p+ = ^{-ttikEs + itikEd - q^ + m\] 



(3.6) 



q -p- = -{-mKEs -ttikEd - q +rnj^^ 
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The form factors F^"^ and have the opposite property under interchange of pion 
momenta and consequently they do not interfere in dV/dq^. Neglecting terms in chiral 
expansion of Oijf') and higher the differential decay rate given in (|3.2|) becomes 



. + _ + dVQ(2) dTp(i) dTp(2) 

— e+e ) = ^ + ^ + ^ . (3.7) 



In fig. ^ we have graphed for each of the three terms on the right hand side of (|3.7|) 



-^— = 2yimK-2m^f-^—^ , 3.8 
dy Tki^ dq^ 

where y = \/q^ / {rriK — Sm^r), the total width of the and we have set = 0. 

Integrating the three terms on the rhs of ( p.7|) over the invariant mass interval > 
{dOMeV)"^ (corresponding to y > 0.14) we find that for wl = 

10^ ■ Br{KL ^ n+TT-e+e-;q^ > (SOMeVf) = 3.8 + 0.78 + 3.4 = 8.0 . (3.9) 

The branching fraction over this range of e~^e~ invariant mass is dominated by the region 
of low q^ and for typical values of ty^ it receives comparable contributions from the form 
factors G and F^'^\ However, in the region of high q^ the branching fraction is likely to be 
dominated by the F^^ form factor. For q^ > (80Mey)'^ (corresponding to y > 0.37) and 
Wl = the three terms on the rhs of ( p.7|) contribute 



10^ ■ Br{KL ^ 7v+n~e+e~;q^ > (SOMeVf) = 0.61 + 0.07 + 1.9 = 2.6 . (3.10) 

A summary of our results for the rate can be found in Table 1. We have displayed the 
contribution to the branching ratio (in units of 10~^) from the three form factors G, F*^^^ 
and F*^^) for different values of the minimum lepton pair invariant mass q^^^- Since the 

(2) 

loop contribution to the form factor F^ is small, it will be difficult to extract a unique 
value for wl from dV /dq^ data alone; a two-fold ambiguity in the value of wl will persist. 
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Lower cut q^^-^ 




Sr(10-8)F(i) 


Sr(10-8)^(.) 


(10MeV)2 


8.8 


3.3 


3.6-3AwL + 0.8wl 


(20MeV)2 


5.6 


1.5 


3.5 -3.3wL + 0.8wl 


(30MeV)2 


3.8 


0.8 


3A~3.2wL + 0.8wl 


(40MeV)2 


2.7 


0.5 


3.1 -3.0wL + 0.7wl 


(60MeV)2 


1.3 


0.2 


2.6-2.4wL + 0.6wl 


(80MeV)2 


0.6 


0.07 


1.9 -l.SwL + 0.4wl 


(100MeV)2 


0.3 


0.03 


1.3-1.2wL + 0.3wl 


(120MeV)2 


0.1 


0.01 


0.74- 0.68wL + 0.16wl 


(180MeV)2 


0.00072 


0.0001 


0.027 - 0.025WL + O.OOGwl 



Table 1: Contributions to the Branching Ratio (10 ^) for a range of q"^^ 



4. The Z-penguin and VF-box AmpUtude 

The short distance VF-box and Z-penguin diagrams give the effective Lagrange density 

^SD = g ^ g7/x(l - 75)^^ 67^756 + h.c. . (4.1) 

Here we only keep the part that contains the lepton axial current (the vector current is 
neglected). It is only the axial current that gives rise to observables that are antisymmetric 
under interchange of e"*" and e~ four momenta, /c-(- ^ k-. 

In (|4.1| ) the quantity ^ receives significant contributions from both the top quark and 
charm quark loops and is given by 

^ = + f ' ^4-2) 

where 

4 = If) + ' (4-3) 

is the sum of the contributions of the Z-penguin and VF-box diagrams. It is convenient 
to express the combination of elements of the Cabibbo-Kobayashi-Maskawa matrix that 
enters in ^ in terms of \Vcb\ and the standard coordinates p + it] of the unitarity triangle 

VZVtd/V:,Vud = {p-l + iv)\Vcb\^ ■ (4.4) 
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A value of \Vcb\ — 0.04 is obtained from inclusive B XcCVf, decay and from exclusive 
B D*eVe decay. Although the values of p and rj are not determined by present data, 
they are expected to be of order unity. 

The quantities and have been calculated including perturbative QCD corrections 



at the next to leading logarithmic level [|T2| [|T^ . There is some sensitivity to the values of 



^QCD 5 "^c and rrit but is of order 10~^ and is of order unity. 

The quark-level Lagrange density in eq. (4.3) can be converted into a Lagrange density 
involving the tt, K and r] hadrons using the Noether procedure. Equating the QCD chiral 
currents with those obtained from chiral variations of the effective lagrangian in eq. ( p.4| ) 
leads to 

CsD = fTrid^EE^T) e7^75e + h.c. . (4.5) 



Expanding out E in terms of the meson fields M we find that the Lagrange density 
implies that the short distance contribution to the Kl — > Tv~^Tr~e~^e~ decay amplitude 
from the VF-box and Z-penguin diagrams is 

M(^^) = + r P^;) u{k-)^,i,vik+) . (4.6) 

5. The Asymmetry Acp 

It is the interference of M^^^^ in (^) with M^^^) in (0) that prod uces the asymme- 



try Acp defined in ( p..8|) . For calculation of Acp it is convenient to use the phase space 
variables used by Pais and Treiman for Ku decay (rather than those used for the total 
rate in Section 3). They are: = (/c_|_ + /c-)^; s = (p+ +p_)^; 9^^ the angle formed by the 
7r+ three momentum and the three-momentum in the tt+tt" rest frame; Oi, the angle 
between the e~ three momentum and the Kl three momentum in the e~^e~ rest frame; 
4>, the angle between the normals to the planes defined in the Kl rest frame by the tt+tt" 
pair and the e~^e~ pair. In terms of these variables 

= sign(sm0) , (5.1) 



\{p- X p+) ■ {k- - k+ 
and the asymmetry is 

r.27r 



^CP = ^rrn J 3 ( [ d(p sign(sin(/>)) [ dc^ dc, ds dq^P X Re (m'^sdY j^{i-y)) 
27(27r)6m|^rK^ Vio J J ^ 



(5.2) 
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where c„ 
are 



cosOtt, Cp = cos^e- The other kinematic functions appearing in this expression 

/3=[l-4m^/.]V2 

X = I — — z — I - sq~ 



1/2 



(5.3) 



In order to evaluate the contributing form factors the foUowing scalar products of four 
vectors are required: 



q-p+ 
q-p- 
P+-P- 



1 



1 



^ {mj^ - s - - -PX cos 

^ (m|- - s-q'^) + ^(3X cos 9.^ 
1 



(5.4) 



;is-2ml 
1 



— - (3X y sq^ sin 6^ sin 6-,^ sin > 



If the variables s and are not integrated over the complete phase space then it is 
understood that the same is to be done for the Kl width in the denominator of 
eq. (13). 

The form factor G does not enter into Re {M^^^^ M'^^^^^ (a sum over e"*" and e~ spins 
is understood). Integrating our cos 6*6 and (j) we find that 



A 



CP 



dcTt ds dq sin6'7r 




Im{i) {Re{F+) + i?e(F_)) + i?e(0 iIm{F+) - Im{F-)) 

(5.5) 

The integration over cos implies that at leading non-trivial order of chiral perturbation 
theory Im{F+) — Im{F-) Im{F^^) — Im{F^^) reflecting indirect CP violation from e 
and Re{F+) + Re{F-) Re{F^^'^) + Re{F^^^) in eq.(p|) . 

Using ( [4.2|) and (|4.4|) we can write the CP violating asymmetry in terms of the real 
and imaginary parts of the CKM elements 

Acp = A^{{p-l)\V,^\%-i,)-A2ri\V,,\''it , (5.6) 

where Ai arises from indirect CP nonconservation (i.e. — mixing ) and A2 arises 
from direct CP nonconservation. We are only able to predict |Acp| since the sign of gg is 
not known. Our expressions for F^^ and F^^ with wl = give (up to an overall sign) 



Ai = 2.7 X 10" 



A, 



3.9 X 10" 



(5.7) 
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for > (SOMeV)^ and 



Ai = 2.4 X 10-^ , A2 = 8.4 X 10-^ , (5.8) 

for g2 > (sOMeV)^. In Table 2 we present Ai and A2 for a range of values of the 
minimum lepton pair invariant mass, g^im normalized to the branching ratios given in 
Table 1 assuming wl = 0. 



Lower cut q^^^ 


Ai 


A2 


(10MeV)2 


2.0 X 10" 


-2 


2.0 X 10" 


-2 


(20MeV)2 


2.5 X 10" 


-2 


3.0 X 10" 


-2 


(30MeV)2 


2.7 X 10" 


-2 


3.9 X 10" 


-2 


(40MeV)2 


2.8 X 10" 


-2 


4.8 X 10" 


-2 


(60MeV)2 


2.7 X 10" 


-2 


6.8 X 10" 


-2 


(80MeV)2 


2.4 X 10" 


-2 


8.4 X 10" 


-2 


(lOOMeV)^ 


2.1 X 10" 


-2 


9.8 X 10" 


-2 


(120MeV)2 


1.8 X 10" 


-2 


0.11 




(180MeV)2 


1.3 X 10" 


-2 


0.13 





Table 2: The CP violating quantities Ai, A2 with tu^, = for different values of q^^ 

We find that direct and indirect sources of CP nonconservation give comparable contribu- 
tions to Acp- In our computation we have neglected final state tttt interactions which are 
formally higher order in chiral perturbation theory. With the values of Ai and A2 given 
in Table 2, \Acp\ is only of order 10~^ and further refinements of our calculation do not 
seem warranted. 

6. The Asymmetry Bcp 

Using the kinematic variables introduced in the previous section the CP violating 
observable Bcp is defined as 

-Bcp =< sign (sin cos 0) > . (6.1) 
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At leading order in chiral perturbation theory it arises from the interference of with 
G^^\ The CP violating form factors F^^^ are not small because they occur at a lower order 
in chiral perturbation theory than the other form factors, F^"^ and G^^\ Consequently, 
as was noted in refs. and 0, Bcp is quite large. Neglecting M^^^^ we find after 
integrating over cj) and cos9e that 

Bcp = , orfo^tff\ r I ^^^^ - ^-)] • (6-2) 

If the variables s and are not integrated over the entire phase space then it is understood 
that the same is to be done to the width in the denominator of ( |6.2| ). The form 
factor G is real at leading order in chiral perturbation theory and the imaginary part arises 
from the phase in — F_ induced by — mixing. The integration over cos 6-,^ implies 
that Fj^ — F_ — F^^ — F^^ in eq. (|6.2| ) . Using our expressions for F^^ and the value of 
1^(2)1 we find that with wl = \Bcp\ ^ 6.3% for > (30 MeVf and \Bcp\ ^ 2.4% 
for q^ > (80 MeV)'^. The asymmetry for a range of values of q^i^ are shown in Table 3. 



Lower cut q"^-^^^ 


\BcP ■ BrilO-^)\ (%) 


(10MeV)2 


134 


(20MeV)2 


78 


(30MeV)2 


50 


(40MeV)2 


33 


(60MeV)2 


14 


(80MeV)2 


6.3 


(100MeV)2 


2.5 


(120MeV)2 


0.92 


(180MeV)2 


0.0086 



Table 3: The CP violating observable \Bcp ■ Br{10 ^)| for a range of values of 

Note that in Table 3 i?r(10~^) denotes the — ^ n~^Tc~e~^e~ branching ratio in units of 
10~^ with the same cut on q^ imposed. We have neglected final state tttc interactions be- 
cause they arise at higher order in chiral perturbation theory. Our prediction for I-BcpI has 
considerable uncertainty because of the neglect of final state tttt interactions and because 
neglected 0{p^) contributions to G seem to be important. 
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7. Conclusions 

In this paper we have calculated the one-photon contribution to the n^n~e~^e~ 
decay rate. We used chiral perturbation theory to determine the form factors and for e~^e~ 
pairs with high invariant mass { q'^ » Ami) found important new contributions that were 
not included in previous work 00 . The amplitude for Kl 7r+7r~e+e~ depends on the 
undetermined (renormalization scale independent) combination of counterterms wl- We 
found that for = (/c+ + /c_)^ > (30 MeV)'^ the branching ratio for n~^n~e^e~ is 

approximately (8.0 — 3.2wl + x 10~^ and for > (80 MeV)'^ the branching ratio 

is approximately (2.6 — I.Swl + OAwj^) x 10~^. 

One interesting aspect of this decay mode is that the CP even component of the 
Kl state contributes at a lower order in chiral perturbation theory than the CP odd 
component. This enhances CP violating effects in —>■ n~^n~e^e~ decay. For example, 
the CP violating observable 0] Bqp =< sign(sin^cos(^) > , where (j) is the angle 
between the normals to the n^n~ and e^e~ planes, is about 6% for q^ > {30MeV)'^ if 
wl = 0. The CP violating observable Acp =< sign(sin^) > arises from the interference 
of W-box and Z-penguin amplitudes with the one-photon part of the decay amplitude. 
Unfortunately, we find that Acp is of order 10~* and hence most likely unmeasurable. 

Chiral Perturbation theory has been extensively applied to nonleptonic, semileptonic 
and radiative kaon decays. The study of Kl 7r+7T~e+e~ offers an opportunity to 
determine the linear combination of coefficients in the 0{p'^) chiral lagrangian that we call 
Wl and to test the applicability of 0{p'^) chiral perturbation theory for kaon decay. 

Some improvements in our calculations are possible. While a full computation of the 
0{p^) contribution to F± and G arising from two-loop diagrams and new local operators 
does not seem feasible it should be possible to calculate the leading contribution to the 
absorptive parts of G and — F_. Note that the absorptive parts come from both 
nn nn rescattering and because of CP nonconservation from nn nivj. We hope to 
present results for this in a future publication. 
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Figure Captions 

Fig. 1. Feynman diagrams contributing to f2^\ 

Fig. 2. Feynman diagrams contributing to the tt^ charge radius, < > , at leading 
order in chiral perturbation theory. 

Fig. 3. The Feynman diagrams contributing to the amplitude for —>■ 7r+7* at leading 
order in chiral perturbation theory. The solid square denotes a vertex from the 
gauged weak lagrangian in ( p.7|) , the solid circle denotes a vertex from the gauged 
strong lagrangian in ( |2.4| ) . Figures 3(a) involve only weak and electromagnetic 
vertices while Figures 3(b) also has a strong vertex. Figures 3(c) are the contri- 
butions from the kaon and pion charge radii (including both loop graphs and the 
tree-level counterterm) . Figure 3(d) is the contribution of the weak counterterm 
as given by ( |2.18| ) . We have not shown the wavefunction renormalization of the 
tree graphs for the process as the sum of these graphs vanish. 

Fig. 4. Feynman diagrams contributing to the CP conserving amplitude for Tr"*" tt" 7* 

at leading order in chiral perturbation theory. The notation is the same as in fig. ^ 
and we have not shown the wavefunction renormalization of the tree graphs for 
the process as the sum of these graphs vanish. 

Fig. 5. The differential decay spectrum as a function of y the invariant mass of the lepton 
pair normalized to rriK — 2m.j^. The dot-dashed curve is the contribution from 
the dotted curve is the contribution from F^^ with wl = and the dashed 
curve is the contribution from G^'^\ The total differential decay rate for wl = 
is given by the solid curve. 
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Figure 2 
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Figure 3(a) 
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Figure 3(c) 




Figure 3(d) 
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Figure 4(a) 
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(m,m) = 



(K+, Jt^) , ( K+ ri ) 
(71+ K^) , (71+ K°) 

Figure 4(b) 
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(m , m ) = ( K , ) , ( K , ri ) 

Figure 4(c) 
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